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Abstract. Recently, an experiment at GSI Darmstadt has observed oscillating decay rates of 
heavy ions. Several controversial attempts have been made to explain this effect in terms of 
neutrino mixing. We briefly describe the experimental results, give an overview of the literature, 
and show that the effect cannot be due to neutrino mixing. If the effect survives, it could, 
however, be explained by hypothetical internal excitations of the mother ions (~ 10~^^ eV). 



The accelerator facility at GSI Darmstadt can produce monoisotopic beams of highly 
ionized heavy atoms and store them for extended periods of time in the Experimental Storage 
Ring. There, an experiment has been performed in which electron capture (EC) decays of 
hydrogen-like ^Pr^*^ and ^g^Pm^""*" ions have been studied using time resolved Schottky mass 
spectrometry This technique allows to detect changes of the ions' revolution frequencies 
which occur upon EC decay. For a small number of stored ions (< 3), Schottky mass 
spectrometry allows for a measurement of the individual decay times. Repeated measurements 
have led to the distributions shown in Fig. [TJ On top of the expected exponential behavior, 
there is a superimposed oscillation with a period of T 7 s. The experiment avoids many 
systematical errors because it provides a quasi-continuous monitoring of the ions. Statistical 
fluctuations are excluded as the origin of the effect at the 99% confidence level 

Several authors have attempted to relate the anomaly to neutrino mixing [2, 0, 0, 0], 
but these attempts have been refuted 0, H, 0] . In the following, we will first argue why the 
GSI experiment is for general reasons distinct from a neutrino oscillation experiment, and then 
present a careful treatment of the problem using the density matrix formalism and the wave 
packet technique. We will show that the GSI anomaly cannot originate from neutrino mixing, 
but could be explained by hypothetical internal excitations of the mother ions (~ 10~^^ eV). In 
the end, we will give a critical overview of the existing literature on the subject. 

Let us start with a general argument why the GSI anomaly is distinct from neutrino 
oscillations, a fact which is not always correctly reported in the literature [13] • Firstly, a Feynman 
diagram for neutrino oscillations, e.g. for the V(, to u^. survival channel, has the following form: 




Important are three different ingredients: The production process of the neutrino, which is 
a weak interaction, ensures that the neutrino is produced in a flavour eigenstate {ve in our 
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Figure 1. Rate of EC decays of hydrogen- like ^^'^Pr and Pm ions as a function of the time 
after the injection into the storage ring. Data taken from the plots in [ij]. 

example). The propagation of the neutrino has to be described in terms of mass eigenstates 
because only those have definite momenta (since the mass of a particle is explicitly included in 
its propagator). The flavour eigenstate is a superposition of mass eigenstates, |z^e) = J2i ^e*il^«)' 
where Uei are the corresponding elements of the leptonic mixing matrix. And finally, the 
detection process involves another weak interaction, where the flavour of the neutrino is measured 
again. Oscillations arise because the propagation amplitudes of the three mass eigenstates are 
different, and hence their mixture after some propagation distance is not the same as at the 
time of production. The oscillation amplitude is given by a coherent sum, 

^ee = j;|C/e.|'e-*^^-, (1) 
i 

from which the standard quantum mechanical oscillation formula is easily obtained. 

However, in the GSI experiment, the situation is completely different, since an electron 
neutrino is produced, but there is no second flavour measurement. The neutrino escapes 
undetected, as shown in the following Feynman diagram: 

undetected neutrino 
mass eigenstate 

mother ion ^ ^', 

\ 

daughter ion ^ 

Since Feynman diagrams describe transitions between states of definite energy and momentum, 
the final state neutrino must be a mass eigenstate due to energy-momentum conservation. 
In a hypothetical situation where the energies and momenta of the mother and daughter ions 
are measured with infinite precision, it would even be possible to tell from kinematics which 
neutrino mass eigenstate has been produced. Only this mass eigenstate, i.e. only one Feynman 
diagram, would contribute and the rate would be proportional to |C/eiPi but there would be no 
oscillations. However, this is far from the real kinematical situation in the GSI setup. 

Realistic energy and momentum uncertainties imply that it is not possible to tell which mass 
eigenstate was produced, so that all of them have to be taken into account, and must be treated 
as distinct final states. They contribute to the total rate as an incoherent sum of the sub- 
amplitudes, since each mass eigenstate lives in its own Hilbert space. The total rate is therefore 
proportional to |C^eiP = implying that in principle there cannot be any mixing effects. Of 
course, this cannot change in a quantum mechanical approximation of field theory. 



Let us now discuss this in more detail, using the density matrix formahsm for a proper 
theoretical treatment of the detection process. We have to take into account that the GSI 
detector is sensitive to the daughter ion, but not to the neutrino. In the density matrix formalism, 
the probability for detecting a particular state |n) in an experiment is given by tr{Pp), where p 
is the density matrix of the system, and P is the projection operator onto |n). If one particle 
is not observed, the projection is not onto one particular state, but rather onto a set of states 
{|n, m) : n fixed, m = . . . 00} , where n denotes the degrees of freedom of the observed particles, 
and m stands for those of the unobserved particle. In the GSI experiment, the detection of a 
daughter state |V'Z),A;) is thus described by the operator 
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The sum and the integral run over a complete set of neutrino mass eigenstates {vj) with momenta 
p^. With the density matrix for the time-evolved mother state \iPm), given by p = \ipM){'ipM\, 
the probability for the observation of \'ipD,k) becomes 



Vk = tr 
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(3) 



We see that the sum over neutrino states is incoherent. Therefore, if Vk contains oscillatory 
interference terms, they cannot be due to neutrino mixing, and would occur also in a hypothetical 
model with only one neutrino flavour. All attempts to explain the GSI anomaly in terms of 
neutrino mixing are thus foredoomed. 

It is, however, imaginable that different components of the mother wave packet \iPm) acquire 
relative phase differences during propagation. If several such components could decay into the 
same daughter state \^pD,k]l^j,Pu), they would induce interference terms in Vk- To see under 
which conditions such a mechanism could explain the GSI anomaly, we will now compute the 
matrix element {ijjD.k] i^j, Pi^lV'A/) in the wave packet formalism (ref. [Hl] and references therein). 
We describe the states \iPm) and \iPd) of the mother and daughter ions by Gaussian wave packets 
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where A = M, D. In our notation, poAf and poD are the central momenta of the wave packets, 
and umi o'd are their momentum space widths. Em and Ed are related to pM and pD by the 
relativistic energy-momentum relation. ■0jv/ is defined such that, at the injection time tM, its 
peak is located at xj,/- Similarly, at the detection time to, the peak of V'd is located at X£). The 
coordinate space Feynman rules yield for the amplitude of the decay into \tpD] ^j-,Pu)'- 



i {il^D;i^j,Pu\tpM) 



d^p 



M 



d^p 



D 



(27r)3V2E^ J {2^f^/2EE 



d^xdt 



27r 



EC 



{Em, Ed, Ey) 



■ exp 



• exp 



{PD - Popf 

(pm — Pom) 

Hi 



+ iEd {t - to) - ipD(x 

) 

- - iEm {t -tM)+ «PA/(x - XA/) 



exp \iEyj t — ipy x] 



(5) 



with being the transition amplitude between plane wave states, as computed in 12, [l 

To evaluate this expression, we are first going to compute the pM integral 
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We expand Em = (p'm + "^m)^''^ the exponent up to first order in pM — Pom (tliis 
approximation neglects wave packet spreading, as has been shown in [lH). The prefactor 
l/y/2EM is assumed to vary slowly over the width of the wave packet, and will therefore be 
replaced by its value at poM, which we denote by 1/\/2Eqm- Similarly, we also neglect the 
energy dependence of the matrix element Jvl^*-" . We then have to evaluate 
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• exp [i{pM - PoA/)(x - xm) - ivoA/(t - tu)] , (7) 
where vqm = Pom/Eqm is the group velocity of the wave packet. The result is 

\3/2 2 2 
exp [ - lEoM (t - Im) + ipoAf(x - xa/)] ( ) —?== exp [ - (x - xa/ - vqa/ (t - tu)) ctm] ■ 

(8) 

After evaluation of the pn integral in a completely analogous way, eq. ([5]) becomes 

i{TPD;i^j,pMM) = f^M^Y f d3xdtUejMf^{EoM,EoD,E^) 

■ exp [- iEoM{t - tA-f) + iPOAf(x - XA/) - (x - XA/ - VoAf(t - iA/))^<7M] 

• exp [- iEooit - to) +ipoD(x - x^) - (x - x^ - voD(t - iD))^o-D] exp [iE^^jt- ipi^x]. (9) 

We see that the x-integral, as well as the subsequent t-integral, are Gaussian; a straightforward 
calculation thus leads to 
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Here, we use the notation 
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The group velocities of the wave packets are given by vo_d,om = Pod,oa//^od,oa/- The real factor 
exp[—fj] enforces sufficient overlap of the wave packets, but is non-oscillatory for Gaussian wave 
packets. The complex phase factor exp[i(j)j] is oscillatory, but is irrelevant for the modulus of 
the matrix element appearing in eq. ([3]). 

We will now construct a hypothetical situation in which the GSI oscillations can be explained 
by a quantum mechanical interference effect, namely by quantum beats of the mother ion. This 
possibility has been pointed out previously in Let us assume that the state of the mother 

ion is split into several sublevels I'i/'i/^); and that, for some reason, the production process creates 
the mother ion in a superposition 

|V'M)=5]a„|vS?), (14) 

n 

where the coefficients a„, have to fulfill the normalization condition |anP = 1- With this 
modification, eq. (jlOp turns into 

i {t/jD] iyj,Pu\ipM) oc ^ a„ exp [ - /j''^ + , (15) 
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where /j*^^ and (j)^^^ are defined as in eqs. ()lip and ()12p . but including an upper index (n) 
for the quantities Eqm, Pom, ^om, Ej, p, v, and ym- For simplicity, we have neglected the 
n-dependence of the normalization factors and of the matrix element. Typically, also the wave 
packet overlap factor exp[— /j'^''] will be almost independent of n, so we can safely omit it in the 
following, assuming it to be absorbed in the overall normalization constant. 

Upon squaring | (1/^/5; z^j, pjy|^/^A/) |, we now obtain interference terms proportional to 

exp[i((/)^'^'' — (/'j'"'')]. To simplify these, let us go to the rest frame of the daughter nucleus, 
in which vqd = and poD = 0. Moreover, we will choose ao = gm = a, and we will expand 
^(") ^("^) £j.g^ order in the small quantities 
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Here, is a real parameter that is determined by the details of the production process. If we 
finally neglect terms of C'((j/E'^"q), we find 
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Using the relation x/j — xj\/ ~ vj^™Q (t/j — tjv/), which is a good approximation for sufficiently 
well localized wave packets, the phase factor can equivalently be written as exp[—i{t£) — 
iM)Am^^/2£'^(]]. To explain the GSI oscillations with T ~ 7 s, one would require Am? ~ 
2.2 • 10""^ eV^, which corresponds to \m?^ — m^^^\ ~ 8.4 • 10~^^ eV. As has been pointed out 
in [13], there is no known mechanism that could split up the ground state of the mother ion by 
such a small amount, nor a known reason why the production process should create a coherent 
superposition of the substates. 



Before concluding, let us discuss why our results disagree with those of several other authors, 
who have claimed that the GSI anomaly is a consequence of neutrino mixing. Ivanov, Reda, and 
Kienle [3] perform a calculation in which the amplitudes Aj = ^(^^'^Pr^^'*' — > ^^"Ce^^^ + 
J = 1 ... 3, receive different phase factors. The appearance of these phase factors can be traced 
back to the assumption of a finite domain for the time integral at the Feynman vertex, and to the 
assumption of momentum non-conservation. However, the authors sum the Aj coherently, and 
thus obtain oscillatory interference terms in the decay rate. To match the observed oscillation 
period T ~ 7 s, a value of ^rn^i ~ 2.22(3) • 10~^ eV^ is required for the solar mass squared 
difference, in conflict with KamLAND results. In a later work [l5[, the authors relate this 
discrepancy to loop-induced Coulomb interactions of the neutrino. Moreover, they a pply their 
formalism also to /^"'"-decays [l^. The treatment of the detection process in Refs. 0, 15,[lB| is 
in conflict with our results, which show that the sum over the Aj should be incoherent rather 
than coherent. Similar arguments have been given previously by Giunti by Burkhardt et 

and by Peshkin [9|. (Note that Ivanov et al. have replied to some of Giunti's remarks 
i, [isl].) Moreover, Giunti has shown another problem, namely that the decay rate computed 
does not reduce to the Standard Model result if the neutrino masses are set to zero 
Further explanation attempts for the GSI anomaly are due to Faber 0] and Lipkin 
Both authors employ various kinematical arguments and assumptions, which in both cases yield 
relative phase differences for the Aj. Faber and Lipkin also sum the amplitudes coherently, in 
conflict with refs. 0, H, and with our discussion. 

Finally, Kleinert and Kienle propose an explanation of the GSI anomaly in terms of a 
"neutrino-pulsating vacuum" The authors interpret the Ve emission in EC decay as the 
absorption of a negative energy z/e from the Dirac sea forming the vacuum. They assume these 
negative energy anti-neutrinos to undergo oscillations, and thus come to the conclusion that the 
rate of EC decay should oscillate as well. However, the Dirac sea contains all three neutrino 
flavours, so that, due to unitarity, its charge remains constant over time. This remains true 
even if the local density of z^e neutrino states should be modified by the presence of an atomic 
nucleus, as alluded by Kleinert and Kienle. 

Besides these theoretical works, two more experiments have been performed investigating EC 
decays of stopped ^^^Pm and ^^''Re atoms 1^, 13] (see also [Ht]). None of these experiments has 




found any oscillatory signature, indicating that the anomaly, if physical, must be related to the 
differing features of the GSI setup. 

In conclusion, we have shown that neutrino mixing cannot cause the oscillating electron 
capture decay rate of heavy ions that has been observed at GSI. If the oscillations were confirmed, 
one must think of new physics, and we have presented a possible, though exotic, mechanism 
which explains the anomaly by hypothetical internal excitations of the mother ion (~ 10~^^ eV). 
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